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Template-based searches for gravitational waves are often limited by the computational cost as- 
sociated with searching large parameter spaces. The study of efficient template banks, in the sense 
of using the smallest number of templates, is therefore of great practical interest. The "traditional" 
approach to template-bank construction requires every point in parameter space to be covered by 
at least one template, which rapidly becomes inefficient at higher dimensions. Here we study an 
alternative approach, where any point in parameter space is covered only with a given probability 
77 < 1. We find that by giving up complete coverage in this way, large reductions in the number 
of templates are possible, especially at higher dimensions. The prime examples studied here are 
"random template banks", in which templates are placed randomly with uniform probability over 
the parameter space. In addition to its obvious simplicity, this method turns out to be surprisingly 
efficient. We analyze the statistical properties of such random template banks, and compare their 
efficiency to traditional lattice coverings. We further study "relaxed" lattice coverings (using Z„ 
and lattices), which similarly cover any signal location only with probability 7;. The relaxed A'^ 
lattice is found to yield the most efficient template banks at low dimensions (n < 10), while random 
template banks increasingly outperform any other method at higher dimensions. 



I. INTRODUCTION 

Matched filtering is the optimal linear detection 
method for signals in additive stationary noise. This 
method is widely used in searches for gravitational 
waves (GWs) in the data of ground based detectors 
(LIGO [2 E], GEO H E], VIRGO |], TAMA ^Tj) 
as well as being among those proposed for future space 
based detectors [5]. The majority of signals that are 
being searched for have predictable waveforms but un- 
known waveform parameters. Matched filtering consists 
of processing the data with multiple waveforms (tem- 
plates) each corresponding to a different set of waveform 
parameters. The different templates are spaced based on 
a metric on the space of signal parameters. The metric 
defines a distance measure directly related to the loss in 
the matched filter signal-to-noise ratio for a given tem- 
plate and signal [9] [TUJ HH [H] . Using the metric, a tem- 
plate bank, usually in the form of a lattice, can be placed 
on the space such that the loss between any putative sig- 
nal and at least a single template in the bank is less than 
a predefined maximum value. 

Template placement for gravitational wave data anal- 
ysis has proven to be a complicated and involved pro- 
cedure even in the relatively low-dimensional spaces al- 
ready searched [13l [HI [15l Hg [17l [18] . Here we discuss 
the possibility of adopting a seemingly far less comph- 
cated template placement method whereby we randomly 
position templates within our search space rather than 
placing them on a lattice. 

It has recently been shown [12] that constructing op- 
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timal template banks can be interpreted as an instance 
of the mathematical sphere covering problem, and that 
results from this field of research can usefully be applied 
to template banks. For instance, the hyper-cubic Z„ lat- 
tice covering is known to become extremely inefficient 
at higher dimensions compared to other lattices, in par- 
ticular the A* lattice, which provides a highly efficient 
covering for dimensions up to n < 24. 

In practice, however, constructing lattice-based tem- 
plate banks often turns out to be problematic, due to 
the difficulties associated with adapting lattice cover- 
ings to curved parameter spaces and performing co- 
ordinate transformations to avoid non-constant metric 
components. Furthermore, even the best lattice cov- 
ering becomes increasingly inefficient at higher dimen- 
sions (n > 10, say), which makes this approach increas- 
ingly unsuitable for problems involving high-dimensional 
parameter-spaces . 

A radically different approach to template-bank con- 
struction consists in relaxing the strict requirement of 
complete coverage for a given mismatch, and instead re- 
quire coverage only with a certain confidence < 1. 
This is a natural step for searches that employ statis- 
tical detection techniques which always involve a finite 
false-dismissal probability. Connected to this idea are 
new types of template banks, commonly referred to as 
"stochastic" , which have recently been studied and ap- 
plied by various groups [20l [211 [211 [23 • Stochastic tem- 
plate banks are constructed by randomly placing tem- 
plates on the parameter space, accompanied by a "prun- 
ing" step in which "superfluous" templates, which are 
deemed to lie too close to each other, are removed. 

Here we study an even simpler approach, which we 
refer to as "random template banks" , in order to dis- 
tinguish it from stochastic banks. This method consists 
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of placing the right number of templates randomly, with 
probability density dependent on the metric determinant 
and without any additional pruning steps. Apart from 
the practical advantage of relative simplicity, this allows 
one to analyze the properties of such random template 
banks analytically and in great detail. For example, we 
can explicitly determine the number of templates Nfi 
required to achieve any desired level of coverage confi- 
dence This paper presents the first detailed study of 
the properties of such template banks, and an explicit 
comparison of their efficiency to traditional full-coverage 
lattice template banks. 

Despite their simplicity, random template banks are 
found to achieve astonishing levels of efficiency compared 
to traditional template banks, especially at higher di- 
mensions. They outperform even the highly efficient A* 
lattice in dimensions above n ^ 6 — 7 for covering confi- 
dences rj in the range of 90% — 95%. 

As a by-product of this study, we also analyze the prop- 
erties of "relaxed lattice" coverings, which share a fun- 
damental feature with random template banks: for any 
signal location, the nominal covering mismatch is guar- 
anteed only with probability 77 < 1. This results in a 
coarser lattice and therefore a reduction of the number 
of templates. We find that these relaxed lattices generally 
result in the most efficient template banks at dimensions 
up to n ~ 11, where random template banks start to 
dominate. 

The plan of this paper is as follows: first we review the 
general template-bank problem and traditional lattice- 
based template banks in Sec. |IT] In Sec. |III| we present a 
detailed analysis of random template banks: We calcu- 
late their template densities, and compare them to tra- 
ditional lattice coverings, and we investigate some of the 
relevant statistical properties of random template banks. 
In Sec. II VI we describe a modification to traditional lattice 
template banks, termed relaxed lattice covering. Sec. |V] 
provides a summary and discussion of the results. 



II. "TRADITIONAL" TEMPLATE-BANK 
CONSTRUCTION 

In this section we briefly review some fundamental 
concepts used in constructing template banks, namely 
the parameter-space metric [9l [TOj |24j and lattice cover- 
ings [TTJ [T^ [12] . One key feature of traditional template 
banks is that they require complete coverage of the pa- 
rameter space, i.e. no point in parameter space is allowed 
to be further away from its closest template than a given 
maximal mismatch. 

Consider an n-dimensional parameter space iS„, with 
coordinates {A*}-Lj. Each point A describes a set of pa- 
rameters of a signal model, which we assume to be an 
accurate description of the true signal family s(t; A). As- 
sume we measured data x{t) containing a signal s(t; As) 
in addition to Gaussian additive noise i.e. x{t) — 

n{t)+s(t; As). Typically one constructs a detection statis- 



tic of the data, X{X; x), say, which is a scalar represent- 
ing the probability that a signal with parameters A is 
present in the data x{t). Due to the random noise fluc- 
tuations n(t), X is a random variable, but with the prop- 
erty that its expectation value X{X; Ag) = E[X{X; x)] has 
a maximum at the true location of the signal A = As . We 
can define a notion of mismatch, or squared length of a 
parameter offset AA = A — As, as the relative loss in the 
expected detection statistic due to this offset, i.e. 

m(AA;As) = l-i^^=5.,(As)AA'AA^ + ... , (1) 

where we use automatic summation over repeated indices 
i,j and the metric tensor g^j is defined via Taylor expan- 
sion of the mismatch m in the small offset A A. Using 
this definition of the metric, the proper volume of the 
parameter space iS„ can now be expressed as 

Ks„ = / dV, with dV = y^ d"A . (2) 

where g = det gij is the determinant of the metric gtj . 
A common approach to the problem of parameter-space 
covering is to use a lattice of templates. Template-based 
searches are often computationally expensive due to the 
large number of templates required, therefore much ef- 
fort has gone into identifying the most efficient covering, 
namely the lattice that requires the fewest templates to 
achieve complete coverage of the parameter space [12]. 
Note that constructing lattice template banks in curved 
parameter spaces is highly impractical, and most of the 
following results implicitly assume that the parameter- 
space metric gij is flat, i.e. we can find coordinates in 
which gij is constant (i.e. independent of parameter-space 
location A). 

A parameter-space point A is considered to be "cov- 
ered" by a template A(fe) if its squared distance to the 
template is smaller than the given nominal mismatch to* , 
i.e. 

gij AA(j.) AA^^j < m* , with A\(k) = A - A(fc) . (3) 

This is equivalent to saying that A lies within the n- 
dimensional sphere of radius R = yfml centered on 
the template The construction of efficient (com- 

plete) template banks is therefore an instance of the 
sphere covering problem, which asks for the sphere ar- 
rangement requiring the smallest number of overlapping 
spheres to completely cover an n-dimensional (Euclidean) 
space [T91I25]. 

A key quantity used in assessing the efficiency of a 
given sphere covering is its thickness. The thickness Q 
is defined 25] as the average number of n-dimensional 
spheres (templates) covering any point in the parameter 
space. For a complete coverage, the thickness therefore 
satisfies by definition > 1 (where in practice equal- 
ity can only be attained for n = 1, in higher dimensions 
there will always be some overlap between spheres). For 
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a lattice covering, the thickness can be conveniently ex- 
pressed as 



e = 



Vn m 



n/2 



VA(m*) 



(4) 



where Vn is the volume enclosed by an n-diniensional 
unit-sphere ^, namely 



Vn 



r"/2 



T{n/2+l)'' 



(5) 



and Va is the volume of a fundamental region of the lat- 
tice A, with covering radius R — y/m^. Note that un- 
der a linear rescaling c, lengths change like R' — c R, 
mismatches like m' = m, and lattice volumes like 
= c" Va- Therefore we see from Eq. Q that the 
thickness is a scale-invariant property, characterizing 
the geometric structure of a covering. In particular O is 
independent of mismatch . A special instance of a fun- 
damental lattice region is the Voronoi cell (also known as 
the Wigner-Seitz cell), which is the set of points that are 
closer to a given template than to any other template. 
Let us also define at this point the volume covered by a 
single template as 



Vt 



n/2 



(6) 



In the following it will also be useful to introduce the 
normalized thickness. 



e_ 

Vn. 



(7) 



which corresponds to the number of templates per unit 
volume in the case of m* = 1. Like the thickness Q, this 
is a scale-invariant property of a covering, independent 
of mismatch m*. As shown in [12], the total number of 
templates of a covering can be expressed in terms of 
the normalized thickness as 



N 



(8) 



which shows that the total number of templates is pro- 
portional to the normalized thickness. 

In the following we will focus on two lattices, namely 
the Z„ (hyper-cubic) and the A*^ lattice, known respec- 
tively for their simplicity and covering efficiency. The 
normalized thickness is known analytically for both lat- 
tices, namely 



2" ■ 



Vn+T 



n(n + 2) 
12(n + 1) 



i/2 



(9) 
(10) 



^ We are using the geometers convention with respect to the defini- 
tion of the n-dimensional sphere (or n-sphere) where the 1-sphere 
represents two points on a hne, the 2-sphere is a circle, etc. 



In the following we will mostly use the normalized thick- 
ness for comparing different covering strategies, since it is 
proportional to the total number of templates (Eq. ([S])), 
which is the quantity we wish to minimize in order to 
reduce the computational cost of searching a parameter 
space. 



III. RANDOM TEMPLATE BANKS 

We now investigate the properties of a new type of 
template bank, which we call the "random template 
bank" , which consists of Nji templates placed randomly 
with uniform probability distribution (per proper vol- 
ume) over the parameter space 5„. Note that contrary 
to the lattice template-banks discussed in the previous 
section, nothing in the following requires the parameter- 
space metric <?jj(A) to be constant or flat. The only prac- 
tical assumption we will make for simplicity is that the 
metric curvature radius has to be large compared to the 
covering radius of one template, so we can neglect met- 
ric curvature in the expression for the volume of a single 
template given by Eq. (|6|. Some practical issues arising 
from non-constant metrics will be discussed in Sec. IIIIEI 



A. Number of required random templates N-r, 

Let us select a point As € Sn which we assume to be 
the location of a signal. We assume that the covering 
sphere with radius R ~ ^Jm^ centered on Ag does not 
intersect the boundary of the parameter-space 5„ , which 
allows us to neglect boundary effects in the following dis- 
cussion. Now consider a single randomly placed tem- 
plate with uniform probability distribution (per proper 
volume) . What is the probability that this template does 
not cover ("miss") within m* the point As? This is equiv- 
alent to the probability that the template does not fall 
within the covering sphere centered on As- The proba- 
bility of falling within this volume is Vt/V^^, and so the 
answer is simply 



P(miss(m,)|5„, Afe = l) = l 



Vn m 



n/2 



(11) 



If we were to place Njz templates randomly in this way, 
the probability that none of the templates cover the point 
As is therefore 



P(miss(m*) |5„, Nn) 



1 



Vn m 



n/2 ' 



(12) 



since in this construction each template location is in- 
dependent of all previously placed templates. It follows 
that the probability that the point As is covered ("hit") 
by at least one template is 



P(hit(TO,)|5„, Afe) = l- 1- 



Vn m 



n/2 ' 



(13) 
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This expression shows that the probabiHty of an (un- 
known) signal location As G Sn being covered by this 
random template bank is always < 1 and we would re- 
quire Nji —^ oo templates to achieve certain complete 
coverage. However, we can relax the requirement on 
certainty of coverage and instead ask how many ran- 
domly placed templates do we need in order to obtain 
a probability rj that an (unknown) signal location As 
would be covered. This is simply given by the solution 
to P (hit(TO*) I Sn, Njz) — rj, which yields 



Afe(?7, TO*, Sn) 



ln(l - 77) 



In ( 1 - to:/' Vn/Vs^, 



(14) 



In practice we will mostly be interested in "large" param- 
eter spaces, in the sense that the parameter-space volume 
Vg^ is very large compared to the volume Vr of one tem- 
plate, and we can therefore Taylor-expand Eq. ( 14 1 in 

the small quantity to*'/' Vn/Vs,^ ^ 1, which yields 



Nji{ri,m^,Sn) 



1 



In 



1 



-n/2 ^ , 



(15) 



where the neglected higher-order terms in the expansion 
correspond to corrections to Nji of order C'(ln(l — rj)/2). 
We see that for a given parameter-space volume V^,^ , the 
two parameters to, and rj completely determine the num- 
ber Nji of random templates we need to place randomly 
on the parameter-space Sn- We therefore introduce the 
notation (to*) to denote an n-dimensional random 
template bank with nominal mismatch to* and covering 
confidence 77. An illustrative example of a 2-dimensional 
random template bank ^■^TZ2{0.1) is shown in Fig. [ijfor 
a nominal mismatch of m* = 0.1, covering confidence 
77 — 0.9, and a parameter-space volume V^,^ requiring 
Nji = 100 templates. By comparing Eq. (15 1 to the gen- 
eral expression in Eq. 
normalized thickness 
namely 



we can directly read off the 
^Tzyq) of a random template bank. 



9^(r;) = In f-^ 

Vn 



(16) 



and from Eq. ([t]) we obtain the corresponding thickness 



(17) 



This expression reveals a very special property of random 
template banks compared to any lattice covering (e.g. 
see Eq. ^ and Eq. (lOl), namely the thickness Q-ji only 
depends on the covering confidence ry, and is independent 
of the dimension n. 

Figure [2] shows a comparison between the normalized 
thickness 9 of the hyper-cubic (Z„) lattice covering, the 
v4* lattice covering and random template banks 
with different covering confidences rj. We see that ran- 
dom template banks beat the efficiency of the A* lattice 




FIG. 1: Example realization of a random template bank 
''■^7?.2(0.1) in n = 2 dimensions (using periodic boundary 
conditions), with nominal mismatch m, = 0.1 and cover- 
ing confidence 77 — 0.9. The parameter-space volume Vs^ 
was chosen such that the resulting number of templates is 
Njz = 100. Template locations are indicated by crosses and 
template boundaries (corresponding to the covering radius 
R — ^/nu) are shown as black circles. The covered and un- 
covered volumes are shaded grey and white respectively. 



(which is the best, or close to the best lattice covering 
currently known for dimensions up to 77, < 24 |T5]) at suf- 
ficiently high dimension n, namely in 77 = 10 for ^'^^TZn, 
77, = 7 for '^■^^TZn, and 77 = 6 for °-^°7^„. Even more 
surprisingly, random template banks beat the theoretical 
lower bound ( "Coxeter- Few-Rogers bound" [23]) on the 
thickness of any covering! This is possible because ran- 
dom template banks do not provide a covering in the 
strict sense, as they leave some fraction of parameter 
space uncovered. 



Distribution of mismatches 



Taking the derivative of Eq. ( 13 1 with respect to the 
nominal mismatch 777*, we obtain the probability den- 
sity function (pdf) for signal mismatches to in a random 
template bank of given Nji and V^^ , namely 



pdf(TO|Afe, Sn) 



nNnVnin'^/^' 



Vn to"/2 



(18) 



which describes the probability of finding a template 
within the mismatch interval [to, to -|- dm] of some lo- 
cation As. This expression is somewhat inconvenient, 
however, as it depends on "extensive" parameter-space 
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relative mismatch in = m/m^ 

FIG. 3: Probability distribution Eq. ( |23[ ) of relative mis- 
matches m = m/mt in random template banks '''^"'TZn of 
dimension n = 1, 2, 3, 6, 12. 
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FIG. 2: Normalized thickness as a function of dimension 
n for hyper-cubic (Z„) and A* lattice covering, and for ran- 
dom template banks ^TZn with different choices of covering 
confidence rj = 0.99, 0.95, 0.90. Also plotted is the Coxeter- 
Few-Rogers bound (CFR), which is a theoretical lower bound 
on the thickness of any strict (77 = 1) covering. 



properties, namely the total number of templates Nfi and 
the parameter-space volume Vs^ . 

In order to rewrite this purely in terms of the "in- 
tensive" parameters to* and 77, let us first note that for 
given nominal mismatch to* and covering confidence 77, 
we can rewrite the required number of templates given 
by Eq. as 



N,z = Qn{v)m-''^^No, with A^o = ^ , (19) 

where we defined Ao as the parameter-space volume V^^ 
measured in units of the unit-sphere volume Vn- Using 
this expression, we can write the probability of a point 
being covered within mismatch to in a random template 
bank ^TZ{m^), using Eq. (13), as 



P(hit(TO) |''7^„(TO*)) ^ 1 



/2 \ AToOtc m. 



1/2 



(20) 

and using the assumption of a large parameter-space vol- 
ume Vg^ , namely Ao S> 1 , we can express this as 



P(hit(m) |''7^„(m*)) w 1 - e" 



-On m" 



(21) 



of the nominal mismatch to*, i.e. 



(22) 



Note that we obviously find P(hit(TO*) | ''7?,„(to*)) = rj. 
Contrary to Eq. (20 1, this expression for P(hit) only de- 
pends on the parameters to and 77. The interpretation 
of Eq. (21 1 is: given a random template bank (m*) 
constructed for nominal mismatch m* and covering con- 
fidence ?7, what is the probability of a signal location 
being covered with a mismatch of at most m. By dif- 
ferentiating this with respect to the relative mismatch to 
we obtain the probability density function of a template 
falling within the mismatch interval [fh, rh + dm] of a 
signal, namely 



pdf(?^|"7^„) = — p(hit(TO) |"7^„(TO*)) 

am 

= |e7jTO"/2-ig-e^rn"/=_ 



(23) 



where we defined to as the mismatch to measured in units 



Figure [3] shows a plot of the mismatch pdf for random 
template banks °-^^7?,„ in different dimensions n. We 
can see in this plot that at higher dimensions the bulk of 
the mismatch probability shifts to the higher end of rel- 
ative mismatches 771. This shows that random template 
banks become more and more efficient as the number of 
dimensions grows. The reason for this is that a bank 
that yields mismatch distributions peaked at lower val- 
ues of the relative mismatch must be more densely pop- 
ulated than a bank with mismatch distribution peaked 
at higher values of the relative mismatch. In order to see 
this it is useful to stress that we construct template banks 
with the only constraint that the maximum mismatch be 
lower than a given fixed value, to* (with probability rj): 
if that constraint is satisfied the bank is acceptable. We 
can construct different banks that satisfy this constraint. 
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with different template densities. In a bank that is very 
densely populated the average distance to the nearest 
template is smaller than in a bank with lower template 
density. Hence in a more densely populated bank the dis- 
tribution of mismatches peaks at a lower value of the rela- 
tive mismatch compared to a lower-density bank. Among 
the two banks, the lower-density one has fewer templates, 
hence it is more efficient. 



C. Spatial parameter-space coverage 

In our discussion of random template banks so far we 
have focused on the probability 77 that an (unknown) sig- 
nal location Ag is covered by a template. If we were 
to construct a number of random template banks, this 
"confidence" therefore describes the fraction of template 
banks in which Ag would be covered. A somewhat related, 
yet different, question is: given a realization of a random 
template bank ^TZnim^), what is the fraction C of param- 
eter space that is actually covered? This spatial coverage 
C, as a property of an individual random-template-bank 
realization, is a random variable, and in the following we 
will analyze some of its relevant statistical properties. 

For a given random-template-bank realization, we de- 
fine the function /(A) on the parameter 



/(A) = 



1 if A is covered 
otherwise. 



(24) 



describing whether a point A is covered ("hit") or not 
( "miss" ) within mismatch to* . Using this we can express 
the spatial coverage fraction C as 



C 



f{\)dV., 



(25) 



where dV — ^^d^X is the volume element associated 
with d"A. Given any point A G 5„ the expectation value 
of /(A) over an ensemble of template banks is given by 

1 

£^[/(A)] = E^'^(/W = J') 

= P(hit(m,) |''7^„(TO*)) 

= V, (26) 

where we used the fact that the probability of any point 
A being covered in ''7?.„(m*) is by construction given by 
the covering confidence rj. Using this we can express the 
expectation value of the spatial coverage C as 

E[C]^-^ [ E[f{X)]dV^rj, (27) 

showing that the expected spatial coverage C is equal to 
the covering confidence 77. Note that this result holds 
true despite the obvious existence of correlations in /(A) 
between neighboring points. 



However, the question of the variance in spatial cov- 
erage C over an ensemble of random template banks is 
substantially complicated by these spatial correlations, 
as we need to evaluate the following integral 



Var [C] 



1 



dv dv'E[f{x)f{y)]\-7r 



(28) 

Since we have no good handle on the spatial correlations 
in /(A) for a given realization of ^TZn, we try to find 
some reasonable approximations. First we approximate 
the above integral as a discrete sum over finite volume el- 
ements AV. In addition we assume that /(A) in each vol- 
ume element AV is completely uncorrelated with /(A') 
in any other volume element, and that /(A) is perfectly 
correlated within each volume element. Applying these 
approximations yields the following estimate for the vari- 
ance: 



Var [C] 



AV 



77(1 - 7?) . 



(29) 



We have found that a reasonably good semi-empirical 
"guess" for the number A^i of independent "uncorrelated" 
volume elements in 5„ seems to be 



Nir^2nNji, and AV 



(30) 



Note that we do currently not have a good theoretical 
understanding of this number, however, it seems to yield 
reasonably good quantitative agreement with the Monte- 
Carlo results on the coverage, as well as the worst-case 
mismatch discussed in the next section. We therefore 
present the analytic estimate as a useful indicator of the 
general behavior and trends at higher dimensions. Sub- 
stituting this into Eq. ( 29 1 , we obtain the following rough 



estimate for the coverage variance 



Var[C] 



77(1-77) 



2nN,. 



(31) 



A key feature to notice is the inverse proportionality of 
the variance to the number of templates Nji and the 
parameter-space dimension n. In Fig. [4] we show example 
distributions of the spatial coverage C, obtained through 
numerical simulations, and compare them to Gaussian 
distributions of mean rj and variance given by Eq. (31 1. 



Each simulation used 2000 random-template-bank real- 
izations (parameter space volumes were computed as a 
function of Nji and 77 and edge effects removed using pe- 
riodic boundary conditions). A Monte-Carlo integration 
using 5 X 10''' points was then performed on each real- 
ization to compute the spatial coverage C. These simu- 
lations become very computationally intensive as the di- 
mensionality of the space increases and hence we limited 
our test scenarios to n < 6 and Nji = 10^. The Gaussian 
distributions using mean E[C] and variance Var[C] seem 
to agree quite well with the results from the Monte-Carlo 
simulations for the range of dimensions considered. 
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FIG. 4: Results of Monte-Carlo simulations of spatial cov- 
erage C for realizations of °'^TZn random template banks in 
n — 1, . . . , 6 dimensions, using Ntz = 10* templates. The his- 
tograms show the distribution of measured coverages C, the 
solid black curves show a Gaussian distribution of mean rj 
and variance given by Eq. ( |31[ ), and the dashed black lines 
indicate the mean values of the measured distributions. 



D. Expected worst-case mismatch 

A related question to the spatial coverage C is the 
worst-case mismatch found in an individual realiza- 
tion of (m*), i.e. the "deepest hole" in the template 
bank. As seen in the previous section, a fraction 1 — C 
of the parameter space iS„ will have mismatch m > m^, . 
Within this uncovered fraction there will be a worst-case 
location with the largest value mw of mismatch. Similar 
to C, the worst-case mismatch TOw is a random variable 
dependent on a given realization of (m*), and in this 
section we analyze the statistical properties of over 
an ensemble of random-template-bank realizations. 

The geometrical definition of the minimal mismatch 
m(A) in any parameter-space point A for a given random 
template bank can be expressed as 

"^(A) = mm [g,, AA^(,) AA^,, | , (32) 

where k indexes the Nji random templates Ajj.) and 
^-^(fc) = ^ A(fe). The worst-case mismatch TOw on the 
space Sn is therefore 

(33) 



mw = max{m(A)} 



Note that this is in fact the definition of the "covering 
mismatch" (or covering radius, squared) of a set of tem- 
plates 25]> ^iid for lattices and other complete (77 = 1) 
coverings, this will be equal to the nominal mismatch, 
i.e. = rrisf. In the case of random template banks, 
however, we control the nominal mismatch to*, but the 



worst-case mismatch could take on arbitrarily large val- 
ues Tow > . One very interesting approach to answer- 
ing this question would be to use the fact that all "local" 
worst-case mismatches occur at Voronoi cell vertices 
(being the so-called "holes" of the covering). For a given 
number of random templates (or "Voronoi seeds" ) the ex- 
pectation value of the number of vertices is known [57^. 
This would be a good start to defining a finite set of 
points for which to draw random mismatch values. The 
problem with this approach, however, is that the vertex 
density grows exponentially with dimension, and vertex 
mismatch values therefore will become highly correlated 
due to their relative "closeness" . In addition, while we 
know the mismatch distribution (see Eq. ( p3| ) at ran- 
domly selected points A, the mismatch distribution on 
Voronoi vertices (being very special points) is different, 
and a complicated function of pdf(TO|''7?,„). 

We therefore employ a more direct and crude approach, 
namely using the semi-empirical "guess" Eq. ( 30 1 for 



the number A'l of statistically "independent" locations 
in Sn, together with the known mismatch distribution of 
Eq. ( 23 1 . A linear rescaling of the whole template space 
will affect all mismatches equally, and therefore it will 
be more useful in the following to consider the relative 
worst-case mismatch TOw = m^jm^, i.e. TOw measured in 
units of the nominal mismatch to* . 

The probability that the largest (relative) mismatch 
TO of Ni independent trials falls within the interval 
[tow , TO-w + dfn] can be expressed as 



pdf(TOw|?7, Ni) dm ■ 



(34) 



where pi is the probability of the mismatch in a single 
trial falling within [tow, Ww + dm], i.e. 



pi = pdf(TOw| "^Tln) dm , 



(35) 



and Pq is the probability of the mismatch in a single trial 
falling within [0, TOw], i-e. 



Po = p(iiit(TOw)r7e„(TO*)) 



(36) 



Using Eq. ( 23 1 this can be rewritten as 

d 



pdf(TO„r'7^„, Ni) 



dm 



P( hit I TO, Tj)^^ 



euN,mZ/'-'e-^-"^"-'' ( 



1 - e 



Afr-l 



(37) 



Using this expression we can easily compute quantiles of 
the worst-case mismatch distribution, namely from 



P{m„ < m\ri, Ni) = / pdf (to^I??, A^i) ^to^ 
Jo 

= fl~e®«""^')^\ (38) 
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FIG. 5: Results of Monte-Carlo simulations of (relative) 
worst-case mismatch mw for realizations of '^'^TZn random 
template banks in n = 1, ... ,6 dimensions, using Nn = lO** 
templates. The histograms show the distribution of measured 
(relative) worst-case mismatches mw, the solid black curves 
show the estimated distribution Eq. ( |37[ ) and the dashed black 
lines indicate the mean values of the measured distributions. 



SO for example we obtain the median worst-case mis- 
match as 



•50% 



1 



e 



In ( 1 - 2-1/^' 



2/n 



(39) 



Using the semi-empirical guess, given in Eq. (30 1, for 
the number A^i of "independent" parameter-space points 
in iS„, we obtain quantitative estimates for the distri- 
bution of the worst-case mismatch. Figure [5] shows a 
few example worst-case mismatch distributions, compar- 



ing the analytical estimate of Eq. (37) to the results of 
Monte-Carlo simulations. These simulations were run 
with 2000 random-template-bank realizations of a ^ '^TZn 
bank containing Nn = lO'* templates each. The pre- 
dicted distributions differ slightly from the simulations, 
but this should be expected given our crude estimation 
of A^i. With the exception of n = 1, the analytical ex- 
pression tends to slightly underestimate the worst-case 
mismatch values. Nevertheless, our rough estimate seems 
to capture the overall trend to smaller values of E[fh^] 
and Var[TOw] with increasing dimension n. This behav- 
ior is also illustrated in Fig. [6) where we have plotted 
the analytic expectation value and variance (computed 
from Eq. (371) of function of dimension n for 



various covering confidences 77 and numbers of templates 
Njz. We see that in higher dimensions the expectation 
value of the relative worst-case mismatch asymptotically 
approaches unity. We also see that as the number of tem- 
plates Nji increases, corresponding to a larger number of 
mismatch "trials" , we obtain an increase in the expected 
value of m^. In addition, we see that the standard devi- 
ation (t[to„] = •y/Var[TOw] decreases both with increasing 
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FIG. 6: Statistical properties of the worst-case mismatch m„ 
in random template banks as function of dimension n, cov- 
ering confidence rj and number of templates Nn- This plot 
shows the expectation value iJ[mw] and standard deviation 
a[mw] derived from the theoretical estimate Eq. (37 1. The 
upper panel shows the case for Nn = 10* templates, while 
the lower panel is for Nn = 10*. 



number of template Nji as well as with higher dimension 
n, so that worst-case mismatches become both smaller 
and more tightly constrained in these limits. For exam- 
ple, as shown in Fig. [5) for n = 6 our Monte-Carlo sim- 
ulations give the mean value of as 1.81 with a stan- 
dard deviation of 0.06 for a random template bank con- 
taining = 10^ templates with a covering confidence 
?7 = 0.9. In using such a template bank one would there- 
fore expect with > 99% (3cr) confidence that the largest 
mismatch is smaller than 1.99 m*. We should also note 
that despite our rather crude estimate for the number of 
independent parameter space locations A^i our model is 
able to estimate the mean value of to within w 10% 
of the values obtained from our simulations. This is in 
fact one of the key results from this investigation, that 
although random template banks by construction only 
provide incomplete coverage, C < 1, the actual worst-case 
mismatch in the uncovered regions is not expected to 
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be substantially larger than the nominal mismatch to* 
especially at higher dimensions. 



E. Practical issues in curved parameter-spaces 

As noted in the beginning of the section, the ran- 
dom template bank results apply without modification 
to curved parameter spaces and non-constant metrics. 
In order to construct such a template bank in practice, 
however, we need to generate a uniform random sam- 
pling in proper volume, which for non-constant metrics 
can be non-trivial. In order to see how to take account 
of non-constant metric components 5ij(A) in the random 
template placement, we note that Eq. ([Ts]) specifies a 
constant uniform probability density p-jz of templates, 
namely 



Ptz 



Nn _ dNjz _ 1 dNjz 
Vs^ ^ ~dV' ^ 7^ 



(40) 



In other words, the non-constant template "pseudo- 
density" pk(A) in coordinate-space satisfies 

pniX) ^ ^ = v^PK - v^^^m:"/' , (41) 

which specifies the required random sampling density in 
coordinate space. 

There are various sophisticated and efficient methods 
for sampling from non-uniform distributions, Markov- 
Chain-Monte-Carlo (MCMC) methods, importance re- 
sampling, and rejection sampling are a few examples. In 
each of these sampling methods it is sufficient to know 
the density pk(A) (and hence the metric determinant) 
only up to some normalising constant factor. However, 
one must know Nji, the total number of random tem- 
plates to draw, which requires accurate knowledge of the 
proper volume, defined by Eq.[2] and therefore ultimately 
one must also have accurate knowledge of the metric de- 
terminant. 

The simplest method, applicable for slowly varying 

template densities is to decompose the parameter space 

(i) 

Sn into smaller patches Sn , which are small enough so 
they can be approximated by a constant metric, and sam- 



pling it uniformly with template density given by Eq. (41 ) 
evaluated at the centre of each patch. 



IV. RELAXED LATTICES 

In the previous section we saw that random template 
banks ^TZn will outperform any covering at sufficiently 
high parameter-space dimension. One of the key fea- 
tures of random template banks, however, is that they 
do not actually provide a strict covering: any point is 
only covered with probability 77 < 1. This allows random 
template banks to beat even the theoretical (Coxeter- 
Few-Rogers) lower bound on the thickness of coverings. 



In higher dimensions it seems to get extremely expensive 
(in terms of number of templates) to cover the "last few 
percent" of a parameter space. Relaxing the requirement 
of complete coverage therefore allows enormous gains in 
efficiency. 

We can now apply this insight to lattice coverings, by 
relaxing the strict "minimax" prescription and instead re- 
quiring a mismatch to* only with probability 77 < 1 . This 
allows us to use a larger maximal mismatch TOmax > ''n* 
for the n-dimensional covering lattice A„ (TOmax), thereby 
reducing the required number of templates. The relation 
between these quantities is given by 



pdf(TO|A„ (TOmax)) dm, 



(42) 



where pdf (to|A„ (TOmax)) is the probability distribution 
of mismatches m for sampled points within a lattice 
An (TOmax)- A Uniform linear rescaling of the whole 
template space will affect all mismatches equally, and 
therefore it will be more useful to introduce the rela- 
tive mismatch to = TO/TOmax, which is invariant under 
re-scalings. Note that contrary to Sec. |III[ here we use a 
relative mismatch defined with respect to TOmax of the lat- 
tice, while the nominal mismatch to* is a-priori unknown 
and will be determined from Eq. (42). In Sec. Ill a rela- 
tive mismatch to = to/to, was used, because there was no 
strict maximal mismatch in this case and we directly pre- 
scribed the nominal mismatch to*. The probability dis- 
tribution of TO is pdf(TO|A„) = TOmax pdf (TO|A„(TOmax)), 

such that r, pdf (to|A„) dfri = 1. We can therefore re- 
state Eq. ( 42 1 in the more useful form 



pdf (to|A„) dm , 



(43) 



where we defined the "effective" relative mismatch fh^, = 
TO* /TOmax < 1, which is determined by the lattice mis- 
match distribution pdf(TO|A„) and the covering confi- 
dence T]. We define the linear relaxation factor r for the 
relaxed lattice ''A„ as 



rCA, 



^J^max 1/2 -, 

= TO* > 1 . 



(44) 



For a given nominal mismatch to* and covering confi- 
dence T], the relaxation factor determines the maximal 
covering mismatch TOmax of the lattice as 



TOmax = TO* CAn) , 



(45) 



which defines the "relaxed lattice" as ''A(to*) = 
A(TOmax)- The normalized thickness 0\^(rj) of a relaxed 
lattice ''An is thereby reduced to 



M = 



r" C^An) ' 



(46) 



with respect to the thickness = (ry = 1) of a tra- 
ditional covering lattice. Note that contrary to random 
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FIG. 7: Results of Monte-Carlo simulations (each using 10® 
points) of the distribution of relative mismatches fh in hyper- 
cubic (Z„) lattices in dimensions n — 1, 2, 3, 6, 12. 
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FIG. 8: Results of Monte-Carlo simulations (each using 10® 
points) of the distribution of relative mismatches m in A* 
lattices in dimensions n = 1, 2, 3, 6, 12. 



template banks, the spatial coverage fraction C of relaxed 
lattices is not a random variable and corresponds exactly 
to the covering confidence, i.e. C = rj and Var[C] — 0. 
Furthermore, the worst-case mismatch is also exactly 
known, namely = rrij^ax- 



As seen in Eq. (43 1 and Eq. (44 1, the relaxation fac- 
tor rCAn) is determined from the mismatch distribu- 
tion pdf(m|A„). The distribution of mismatches depends 
strongly on the type of lattice A„ and the dimension n. 
Unfortunately, these mismatch distributions are gener- 
ally not known analytically, and we need to resort to 
Monte-Carlo simulations to determine them. 

In order to sample the lattice mismatch distribution, 
we uniformly pick points in parameter space, then find 
the closest lattice template and determine its relative 
mismatch to — m/m^i^^. For finding the closest tem- 
plate, we use an elegant and efficient method described 
in [28j, which is based on the "fast quantizing" algo- 
rithms available for many well-known lattices (see [29] 
and [2S^ Chapter 20). Here we focus on two lattices only: 
the simple but inefficient hyper-cubic lattice Z„, and the 
highly efficient covering lattice A* . Some resulting sam- 
pled mismatch distributions are shown in Figs. [7] and [8] 
for different dimensions n. The Monte-Carlo simulations 
used 10^ sampling points for each lattice A„ (except for 
A* in n = 18, 19, where 10^ points were used), and the 
mismatches were binned into 1000 mismatch bins. The 
"jitter" seen in Figs.[7]and|8]illustrates the intrinsic sam- 
pling ffuctuation in these simulations. 

The resulting relaxation factors r CAn) for covering 
confidences -q = 0.99, 0.95, 0.90 obtained via Eq. ( [44 ) 
are shown in Fig. [9] The errors on the relaxation factors 
were determined using a Jackknife estimator (see [55], 
using 100 subsets) and are found to be below 0.04% in all 
cases. We see in Fig. [9] that the hypercubic (Z„) lattice 
can be relaxed substantially more than the ^* lattice, 
which is also apparent from the pdfs in Fig. [7] and [8] the 
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FIG. 9: Linear relaxation factors K) for Z„ and AJj lattices 
and covering confidence r) — 0.90, 0.95, 0.99. 



mismatch distribution of Z„ is much more "wasteful" , 
as it increasingly concentrates around fh = 1/3. The 
A* lattice on the other hand, which is a highly efficient 
covering lattice, has the bulk of mismatches concentrated 
closer to the maximal mismatch m = 1. 

In Fig. [To] and Table [T| we show the resulting covering 
thickness 6'a„ (^7) of the relaxed and ''A* lattices in 
comparison to random template banks ''7?.„. We see that 
while relaxed lattices are substantially more efficient than 
traditional complete-coverage lattices, at higher dimen- 
sions the random template banks eventually still outper- 
form them. At low dimensions, n < 10 say, the relaxed 
''A* lattice provides the most efficient covering we have 
found so far. However, having studied only relaxed Z„ 
and A* lattices so far, it is conceivable that other lattices, 
while not necessarily very good covering lattices, could 
provide an even better relaxed lattice than ''A*. More 
work is required to study this possibility. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 




0.50 


0.50 


0.65 


1.0 


1.7 


3.4 


7.1 


16 


38 


98 


261 


729 


2e3 


6e3 


2e4 


7e4 


2e5 


8c5 


3e6 


1 n i * 
l.U ^ * 


0.50 


0.38 


0.35 


0.36 


0.40 


0.49 


0.65 


0.90 


1.3 


2.1 


3.3 


5.6 


9.9 


18 


34 


65 


130 


266 


559 




0.49 


0.43 


0.43 


0.49 


0.63 


0.86 


1.3 


2.0 


3.2 


5.5 


10. 


19 


36 


73 


149 


322 


699 


2e3 


4e3 


0.99 A* 


0.49 


0.35 


0.29 


0.28 


0.30 


0.35 


0.44 


0.59 


0.83 


1.3 


2.0 


3.2 


5.5 


9.7 


18 


33 


65 


131 


268 


99^7"? 


2.3 


1.5 


1.1 


0.93 


0.87 


0.89 


0.97 


1.1 


1.4 


1.8 


2.4 


3.4 


5.1 


7.7 


12 


20 


33 


56 


99 




0.48 


0.36 


0.33 


0.35 


0.42 


0.54 


0.75 


1.1 


1.8 


2.9 


5.0 


8.9 


16 


32 


63 


130 


274 


590 


le3 


0.95 A-* 


0.48 


0.31 


0.26 


0.24 


0.25 


0.29 


0.36 


0.47 


0.66 


0.98 


1.5 


2.5 


4.2 


7.3 


13 


25 


47 


95 


192 


0.95'!-) 


1.5 


0.95 


0.72 


0.61 


0.57 


0.58 


0.63 


0.74 


0.91 


1.2 


1.6 


2.2 


3.3 


5.0 


7.9 


13 


21 


36 


64 




0.45 


0.31 


0.28 


0.28 


0.33 


0.41 


0.55 


0.80 


1.2 


2.0 


3.3 


5.8 


10 


20 


38 


77 


160 


339 


748 


0.90 4* 


0.45 


0.29 


0.23 


0.21 


0.22 


0.25 


0.31 


0.40 


0.56 


0.82 


1.3 


2.0 


3.4 


5.9 


11 


20 


38 


75 


150 




1.2 


0.73 


0.55 


0.47 


0.44 


0.45 


0.49 


0.57 


0.70 


0.90 


1.2 


1.7 


2.5 


3.8 


6.0 


9.8 


16 


28 


49 



TABLE I: Normalized thickness 6 in dimensions n < 19 for traditional lattice covering (^"Z„, ^'^A^), relaxed lattice covering 
CZn, ''A*) and random template banks (^TZn) with covering confidences rj — 0.99, 0.95, 0.90 respectively. Boldface indicates 
the lowest thickness at given covering confidence rj and dimension n. 
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FIG. 10: Normalized thickness 9 as function of dimension n 
for strict lattice covering A„, relaxed lattices ''A„ and random 
template banks ^TZ„ for values = 0.99, 0.95, 0.90 of covering 
confidence. The upper panel is for the Z„ lattice, while the 
lower panel shows the case of the A* lattice. 



V. DISCUSSION 

Our results show that giving up deterministic certainty 
of coverage of a parameter space can result in large gains 
in efficiency, by substantially reducing the number of 
required templates. The prime example of such a re- 
laxed covering is the random template bank construction, 
which we defined as templates placed randomly with uni- 
form probability distribution (per proper volume) over 
the parameter space. Such a random template bank 
''7?,„(m*) covers any signal location (excluding bound- 
ary effects) within mismatch with probability r/ < 1. 
We have found that the template density of these ran- 
dom template banks can be significantly lower than that 
of even the most efficient (complete) covering, and that 
this advantage increases in higher parameter-space di- 
mensions. 

The exclusion of boundary effects in the random tem- 
plate bank analysis is valid for situations in which the 
projected length scale of the templates in each of the 
search dimensions is ^ the width of the parameter space 
in the corresponding dimension. We should note that 
boundary effects have also been excluded in the calcula- 
tion for the lattice coverings (traditional and relaxed) and 
in practice boundary effects will be equally problematic 
for all template bank strategies. 

Other studies have recently started to investigate a 
somewhat different random template placement strategy, 
referred to as "stochastic" template banks [201 HIl [HI H3] ■ 
The key difference of these methods is that they involve 
a "pruning stage" in the random template placement, 
which is aimed to remove templates that are too close to 
each other. Albeit not yet completely quantified, these 
methods could potentially produce even more efficient 
template banks for equivalent coverage. The disadvan- 
tage, however, is that the template placement procedure 
is more complicated and more computationally intensive. 
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requiring the comparison of each template with every 
other. Stochastic template banks might therefore be less 
suitable for higher dimensions and larger numbers of tem- 
plates. In addition, it seems unclear how much efficiency 
can be gained by such a pruning step, especially in higher 
dimensions. 

Relaxed coverings do not provide complete coverage of 
the template parameter space. Applying such a scheme 
therefore affects the final results as an additional uncer- 
tainty. For gravitational-wave searches such uncertainty 
can be well maintained at a level comparable to other 
uncertainties of the problem - typically of the order of a 
few percent - hence not significantly affecting the overall 
degree of confidence of the result. 

In computationally limited searches, relaxed template 
banks allow significant reductions in computational cost. 
Through the re-investment of this saved computational 
cost, this can yield an increase in sensitivity and breadth 
of the search. 

We have investigated some of the relevant statisti- 
cal properties of random template banks, in particu- 
lar the spatial parameter-space coverage fraction C and 
the worst-case mismatch rriw expected in individual 
template-bank realizations. We have performed Monte- 
Carlo simulations to determine the statistical distribu- 
tion of these quantities, and we have found a rough ana- 
lytical estimate, which shows reasonably good agreement 
with the numerical results. These results show that the 
variance of the spatial coverage fraction C is inversely 
proportional to both the dimensionality of the space and 
the number of random templates. More importantly, the 
worst-case mismatch > m* is found to be typically of 
the same order of magnitude as the prescribed nominal 
mismatch to*, and is rapidly approaching to* with in- 
creasing dimension n. At n = 4, for example, the largest 
mismatch for Nji = 10^ templates and rj — 0.9 has an 
expected value of ~ 3 to* with a standard deviation of 
0.09 TO* (see Fig.[6|. At n = 12, the distribution peaks at 
even lower values, with an expectation value of ^ 1.5 to* 



and a low standard deviation of 0.014 to*. Even though 
the coverage at the nominal mismatch to* only holds in a 
statistical sense defined by the confidence 77, in practice, 
at dimensions greater than n = 11, the worst-case loss 
is within a factor of 1.5 of the nominal value, with very 
high confidence. 

Inspired by these results, we have also investigated the 
properties of "relaxed" lattice coverings, which follow the 
analogous prescription of a nominal covering mismatch 
TO*, achieved with probability 77 < 1 for any signal lo- 
cation. This leads to lattices with larger maximal mis- 
match TOiuax > '71*, thereby reducing the required num- 
ber of templates. We have analyzed the properties of 
relaxed Z„ and A* lattices, and we have found that the 
relaxed ^A'^ lattice provides the most efficient covering 
found so far for n < 10, for covering confidences ry > 0.90. 
In higher dimensions, however, random template banks 
outperform any other method considered so far, includ- 
ing relaxed lattice coverings. 

Possibly the greatest advantage of random template 
banks for "real- world" applications, however, is their 
practical simplicity. Constructing lattice template banks 
is notoriously difficult, especially regarding the handling 
of curved parameter spaces and non-constant metrics. 
Random template banks, on the other hand, are nearly 
trivial to construct, even in spaces with non-constant 
metrics: one only needs to adjust the spatial probabil- 
ity density according to the metric determinant. 
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